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An analytical formulation for the second virial coefficient is given for a spherically symmetric 
potential function of a simple form which has been called a realistic pair potential. This potential, given 
as an inverse relationship with r as a function of U, is capable of being much softer than a Lennard- 
Jones in the extreme of close approach. Differing forms of the result are given, including an expression 
by means of the generalized hypergeometric series. The result is also expressed in terms of the second 
virial coefficient for the (12, 6) potential. 
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One of the uses of a pair potential function is for the 
representation or prediction of the second virial coeffi- 
cient contribution to the equation of state due to the 
two molecules having the specified interaction poten- 
tial energy. Various potential functions have been used 
with some success, the most notable of which may 
be the Lennard-Jones (especially the 12, 6), the 
Morse function, and the exp-6. 

In spite of the varied nature of the different func- 
tions, it is sometimes pointed out that there is so much 
A similarity in second virial behavior for them that it 
is difficult to select one in preference to another on 
the basis of second virials alone. Even if this be true, 
there still may be some interest in yet another poten- 
tial for which the second virial may be formulated 
explicitly, particularly since this formulation can be 
made in terms of already tabulated second virial co- 
' efficients for the (12, 6) function. 

A potential function has been introduced earlier 
[l] 1 in which the value of r/r m is given as a function 
of energy for the usual 1/r 6 type of London interaction 
according to the relation 

(r/r m ) 3 = (/>- 1/2 exp [-3e 1 (c/)-l)-3e 2 ((/>-l) 2 . . .]. 

? 

(1) 

The quantity </> is related to the energy U according to 
the equation 



The range of values of </> is from zero for very large r 
where U is zero, through <j) = l for r = r m where U has 
the value — €, and to very large values where U be- 
comes very large as r becomes very small, according 
to the expected pattern of behavior. 

The second viral coefficient may for the present 
simple case of a spherical potential function be 
written as 



B = N [°° (l-exp(-U/kT)). 27rrWr 



(3) 



uie=(<t>-iy-\. 



(2) 



if the pair is composed of like molecules. Integration by 
parts transforms the integral to 

B= (2tt/3) (NolkT) f * [H>] exp (-U/kT) dU, (4) 

where the quantity [r 3 ] is a function of U and represents 
the simple r 3 value in the positive energy region but 
represents the difference of i 3 for the repulsion side 
of the potential valley minus r 3 for the attractive side 
if the potential energy is negative. This double-valued 
character is removed if the variable (f> is introduced, 
since the integral then is 

B=(27rl3)N (2elkT) 

r* exp [-(<t) 2 -2<t>)elkT] (0- l)d<f> (5) 



/; 



- 
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with (f> advancing simply from zero to infinity and with 
r 3 having its ordinary algebraic meaning. 

If we now introduce the present potential in the form 
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(r/r m ) 3 = 4>-" 2 exp(-a(/> 2 -6</>+c) 



(6) 



where a = 3e 2 , /> = 3ei— 6e 2 , and c=3e r — 3e 2 , we have, 
with 7 1 * for jfcr/e, 

B=(27Tl3)N n rf„(2IT*) f * exp [c- (b-2IT*)<f> 

-(a + l/r*)</) 2 ](</)-l)</)-" 2 ^ (7) 

Expansion of the exponential involving the first power 
of <f> using/3=6 + 2/r* anda=a+l/r* gives 



^-==-|^a- 1 / 4 (7 , **) 3 / 4 ^/2{(l/2 + fa- 1 / 2 7 , **^/ 2 ) 

Dm ' 

B*(T**) + (2 + 4a-" 2 r**" 2 )fif (T**) 
" (12) 



Bib. 



1+aT* 



[2-bT*] 2 \2-bT* 

{[4+(3-2b)T*]B*(T**) 

+ 8[2+ (4-b)T*]Bf(T**) + 16T*BUT**)} (13) 



fl=(2<n-/3)JVo/i(2/r 



*) f * 
Jo 



exp [c — ac/> 2 ] 



Another arrangement of the result can be arrived 
at in the following way. If the second sum in the ex- < 
pression for B/b m is written out in detail through 
(8) several terms, it is possible to recognize that there is 
a relation to the generalized hypergeometric function 
and the confluent hypergeometric function in particu- 
The substitution of x for (j) 2 puts the result in the lar. With r representing (4/r**) 1/2 , the summation 
form 



y £^(4>-i)4>-" 2 ^. 



B= {2TT/3)Nor m (e c IT*) J £| 

rc = U ' 
e -ax( x nl2-\ jAr - x nl2-m)dx 



/; 



(9) 



becomes 

+ r(3/4)[,F 1 (3/4;3/2;T 2 /4)]. 






(14) 



or 



n = 



B= (2irl3)Nor 3 m {e c IT*) 

[ a - n /2-3/4 r(ri /2 + 3/4)-a-' | / 2 - 1 ' 4 r(W2+l/4)]. 



This can be written as#=6 m . (Blb m ) with6,„= (277/3) 
N r m , and 

B l bm= {e c /T * } \ a -3i4 J r(/ l /2 + 3/4)(4/r**)'" 2 //i! 

L n = 



which also has a relationship to Bessel functions 
[2] We further note that if the ordinary hypergeometric 
function F(a, (3; y; f) with €=r)IP is considered and 
the limit taken as /3 _1 -»0, the result is found to be 
n iFi(a; y; 7)). A function equivalent to the original 
'J hypergeometric function according to a standard 
relation is (1 - %) y - a -VF{y — a, y — jS;y;|). A similar 
use of £=r)lf3 and a passing to the limit /3 -1 -*0 
produces an equivalent form e\F t (y — a^; — iq). 
Then for the case with a = iy the" equating of the two 
expressions leads to 



-a" 1/4 £ r(n/2 + l/4)(4/r 



**)«/2/„! 



e-^ 2 i Fi(ir; y; i?) = e^iXiy; y; - 2) ; (15) 



- 



(11) 



where 7 1 ** is an abbreviation for 4a//3 2 . 

Another formulation can be given by which the 
summations may be avoided by various relationships 
according to which the first summation is equal to 

- (V2) (r**) 5 / 4 [(3/8)B*(^**) +4fl*(r**) +2B*(r**)] 

and the second is 

_2_±(7 , **)3/4[B*(7 l **)+4B*(7 , **)] 9 



where 5*, Bf, and# 2 * are the usually tabulated second 
virial coefficients and dimensionless derivatives for 
the Lennard-Jones (12, 6) potential function. 
Using these, (11) can be written 



from which it is evident that both expressions, on 
either side of the equation, are even functions of n. 
It would appear plausible that a series involving even 
powers only would need half as many terms to attain * 
a given quality of approximation as a related one with 
both odd and even powers. Thus, we would presume 
that a more rapidly convergent form for the result 
should exist and might be looked for. For this pur- 
pose we examine the differential equation [2] 



a Y = (y~ x)y' +xy" , 



(16) 



for which y=iFi(a; y; x) is known to be a solution. 
For the special case with a = iy, if we use z=e~ xl2 y 
ory=e xl2 z and substitute into the differential equation, 
we obtain 

xz" + yz f -ixz = (17) 

Solution in series by the method of Frobenius gives 
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: 1 + l-8\4/ + 1-2-8(6+1) W 



l-2-3-8(8+l)(8 + 2) \4/ 



with Y\(t) obtainable by differentiation according to 
Y 1 (T)=d¥ (T)ldr, and with r= (4/7 1 **) 1 / 2 . 

It may be seen that this particular hypergeometric 
U") function has a simple rule for differentiation; namely, 



of, 8; 



x Y 
4 



(19) 



d[ F 1 (a;x)]ldx=a- l i) F l (a+l;x). (22) 

Thus a procedure may be found to evaluate Bjbm 
and any of its derivatives which are of interest. 



where 8=(l+y)/2, according to the notation [3] 
for the generalized hypergeometric series. 

In summary, the final series obtained may be sub- 
stituted into the earlier expressions with the result 
that the second virial coefficient may be obtained from 

B/b m =e c T*- 1 [a-^Y l (T)-a-H*Yo(T)] (20) 



Fo(r) = [r(l/4)oF,(3/4;T 4 /256) 

+ T(3/4)t oFi(5/4; r 4 /256)] exp (r 2 /8) (21) 
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